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ABSTRACT 

The two pulsating pre-main sequence (PMS) stars V 588 Mon and V 589 
Mon were observed by C0R0T for 23.4 days in March 2008 during the Short Run 
SRaOl and in 2004 and 2006 by MOST for a total of -70 days. We present 
their photometric variability up to 1000 fi Hz and down to residual amplitude 
noise levels of 23 and 10 ppm of the C0R0T data for V 588 Mon and V 589 
Mon, respectively. The C0R0T imagette data as well as the two MOST data sets 
allowed for detailed frequency analyses using Period04 and SigSpec. We confirm 
all previously identified frequencies, improve the known pulsation spectra to a 
total of 21 frequencies for V 588 Mon and 37 for V 589 Mon and compare them 
to our PMS model predictions. No model oscillation spectrum with / = 0, 1, 2, 
and 3 p-modes matches all the observed frequencies. When rotation is included 
we find that the rotationally split modes of the slower rotating star, V 589 Mon, 
are addressable via perturbative methods while for the more rapidly rotating 
star, V 588 Mon, they are not and, consequently, will require more sophisticated 
modeling. The high precision of the C0R0T data allowed us to investigate the 
large density of frequencies found in the region from to 300 /iHz. The presence 
of granulation appears to be a more attractive explanation than the excitation 
of high-degree modes. Granulation was modeled with a superposition of white 
noise, a sum of Lorentzian-like functions and a Gaussian. Our analysis clearly 
illustrates the need for a more sophisticated granulation model. 

Subject headings: asteroseismology — techniques: photometric — stars: pre- 
main sequence — stars: variables: delta Scuti — stars: individual (V 588 Mon, 
V 589 Mon) 



1 Based on data from the MOST satellite, a Canadian Space Agency mission, jointly operated by Dynacon 
Inc., the University of Toronto Institute for Aerospace Studies and the University of British Columbia with 
the assistance of the University of Vienna. 

2 The C0R0T space mission, launched on December 27th 2006, has been developed and is operated by 
CNES, with the contribution of Austria, Belgium, Brazil , ESA (RSSD and Science Programme), Germany 
and Spain. 



-3 - 



Introduction 



V 588 Mon (HD 261331, NGC 2264 2) and V 589 Mon (HD 261446, NGC 2264 20) 
are two pre-main-sequence (PMS) pulsating stars for which there exists strong evidence that 
they are members of the young open cluster NGC 2264. Th e proper motions for both stars 
are in agreement with the clusters average proper motion (IHog et al. 1 120001 ) and both fit 
the clusters HR- and color- magnitude d iagrams well. A radial velocity measurement only 
exists for V 589 Mon (IStrom et al. Ill97ll ) but it is consistent with the values for other cluster 
members. Finally, the radial velocities of emission lines in the optical spectra of the two stars 
caused by interstellar gas match the radial velocities of the cl uster suggesting that bo th stars 
are indeed embedded in gas clouds belonging to the cluster (jKallinger et al. Il2008ai ). 



The cluster has a diameter of ~39 arcminutes and belongs to the Mon O B 1 association 
The age of NGC 2264 is reported in the lit erature to lie betw een 3 (e.g., I Walker 1 1 195 6 



Sung et al. 1 120041 ) and 10 million years (e.g., ISagar et al. 1119861 ). With such a young age 



the cluster's main sequence only consists of massive O and B type stars, while stars of later 
spectral types are still in their pre-main sequence (PMS) phase. Therefore, V 588 Mon and 
V 589 Mon having spectral types of A7 and F2, respectively, have not arrived on the zero-age 
main sequence (ZAMS) yet. 



The S Scuti like variability of V 588 Mon and V 589 Mon was first reported by iBreger 

( 119721 ). Hence, they were the first pulsating PMS stars discovered. In the mean time the 
number of known 5 Scuti like PMS pulsators has increased from 36 (IZwintz 1120081 ) to ~ 60 



(Zwintz, K., private communication) due to dedicated observations from ground and from 
space. 

Pulsating PMS stars have intermediate masses, i.e. between ~1.5 and 4 M Q and can 
become vibrationally unstable when they cross the instability region in the Hertzsprung- 
Russell (HR) diagram on their way to the ZAMS. Pre- and post-main sequence evolutionary 
tracks for the same stellar mass intersect several times close to the ZAMS which makes the 
determination of the evolutionary stage of a field star from only its effective temperature, 
luminosity and mass ambiguous. Additional information such as typical observational evi- 
dence for the PMS evolutionary stage (i.e. emission lines, IR excess, an X-ray flux, located 
in an obscured region on the sky etc.) or membership to very young open clusters is needed 
to resolve this ambiguity. Another way to distinguish the evolutionary stages can come 
from the asteroseismic interpretation of the observed pulsation frequencies ( IGuenther et al. 
2007h . 



Time series photometry f or V 588 Mon and V 589 Mon has been obtained from a multi- 
site ground based campaign ( jKallinger et al. I l2008a| ) and from two observing runs of the 
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Cana dian MOST space telescope flWalker et al. 1120031 ) in 2004 and 2006 (IGuenther et al. 



20091 ). The 8 and 12 frequencies common to these three data sets fo r V 588 Mon and V 58 9 



Mon, respectively, were submitted to a first asteroseismic analysis (IGuenther et al. II2009I ). 



The accuracy of the MOST observations is higher than that of the ground based data. Hence 
it is not surprising that the two MOST data sets yield more significant frequencies at lower 
amplitudes that were not found in the ground based observations. The CoRoT observations 
of the two PMS pulsators are new and independent data sets of unprecedented accuracy that 
allow us for the first time to investigate other effects (e.g., granulation) in PMS stars. As 
both MOST data sets and the CoRoT data are available to us, we use them together for a 
detailed pulsational analysis. 



2. Observations and Data Reduction 



2.1. MOST observations 



The MOST (Microvariability and Oscillations of STars) space telescope (IWalker et al. 



20031 ) was launched on June 30, 2003, into a polar sun- synchronous circular orbit with an 
orbital period of ~ 101 minutes (corresponding to an orbital frequency of 14.2 d _1 ). From 
its orbital vantage point, MOST can obtain uninterrupted observations of stars located in 
its Continuous Viewing Zone (CVZ) for up to 8 weeks. The MOST satellite houses a 15-cm 
Rumak-Maksutov telescope feeding a CCD photometer through a single, custom broadband 
optical filter (covering wavelengths from 350 to 750 nm). 

MOST can supply up to three types of photometric data simultaneously for multiple 
targets in its field. The mission was originally intended only for Fabry Imaging, in which an 
image of the entrance pupil of the telescope - illuminated by a bright ta rget star (V < 6) - is 
projected onto the instrument's Science CCD by a Fabry microlens (see lReegen et al. 1120061 ) 
for details. After MOST was operating in orbit, the pointing performance of the satellite was 
improved so much that a new mode of observing, Direct Imaging, was made practical. Direct 
Imaging is much like conventional CCD photometry, in which photometry is obtained from 
defocussed images of stars in the open area of the CCD not covered by the Fabry microlens 
array field stop mask. In the original design, no scientific information was available from 
the Guide Stars used for the ACS (Attitu de Control System), but now precise ph otometry 



is possible for these stars as well (see e.g., [Walker et al. Il2005t lAerts et al. 1 120061 ) . 



V 588 Mon and V 589 Mon were observed in the Direct Imaging mode, first from 
Dec ember 6, 2004, to J anuary 24, 2005, and second within the NGC 2264 observations 



sec 



Zwintz et al. 1120091 ) from December 7, 2006, to January 4, 2007. The light curves have 
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therefore time bases of 48.13 days in 2004 and 22.77 days in 2006. In 2004 on-board exposures 
were 15 seconds long sampled 2 times per minute. A 5-day subset light curve of V 588 Mon 
from 2004 is shown in panel a of Figure [U the respective V 589 Mon 5-day subset light curve 
in panel d. The 2006 data had exposure times of 1.5 seconds, 16 co nsecutive images wer e 



added on top of each other and the sampling time was 24 seconds (IZwintz et al. 1 120091 ) . 



5-day subsets of the light curves for V 588 Mon and V 589 Mon from the 2006 data are given 
in panels b and e of Figure [TJ respectively. 

Data reduction o f the MOST Direct Imaging photometry was conducted using the 
method developed by iRowe et al. I (120061 ) which combines classical aperture photometry 



and point-spread function fitting to the Direct Imaging subrasters. 



2.2. CoRoT observations 



The CoRoT satellite (IBaglin 1120061 ) was launched on December 27, 2006, from Baikonur 



aboard a Soyuz rocket into a polar, inertial circular orbit at an altitude of 896 km. With 
its 27-cm telescope, CoRoT can observe stars within a field of view of about 1.3 x 2.6 deg 2 
located inside two cones of 10-degree radius, one at a right ascension of 06:50 (galactic 
anticenter direction) and the other at 18:50 (galactic center direction). 

The CoRoT space telescope originally had two CCDs devoted to asteroseismology for 
stars with 5.7 < m v < 9.5 mag and two CCDs dedicated to the search for exoplanets 
where ~6000 stars in the magnitude range from 10 to 16 mag in R per CCD are monitored. 
Observations of the young open cluster NGC 2264 were conducted as dedicated Short Run, 
SRaOl, uninterru pted for 23.4 days in March 2008 within the framework of the Additional 



upt 



Programme (AP: I Weiss I \200d) . At the middle of SRaOl the ascending mode of the orbital 
plane was 15.67°. The complete cluster was placed in one of the two Exofield-CCDs and 
data were taken for all stars in the accessible magnitude range. The 100 observed brightest 
stars in the field of the cluster were primary targets to search for stellar pulsations in PMS 
stars. A detailed description will be given in Zwintz et al. (2010, in preparation). 

With V=9.7 mag and V=10.3 mag, V 588 Mon and V 589 Mon are nominally too bright 
to be observed in the CoRoT Exofield. Obtaining data for these two stars was only possible 
using the method of the so-called CoRoT imagettes. 

Per observing run up to 20 imagettes, 15 x 10 pixel large CCD subwindows, can be 
defined in each of the two CoRoT Exo-field CCDs. SRaOl was the first CoRoT observing run, 
for which the imagette data were requested for scientific use and allowed to observe stars that 
normally would be too bright for the Exo-field. Each of the imagettes is submitted to a special 
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data reduction pipeline developed by the Laboratoire d'Astronomie de Marseille (LAM) and 
the Laboratoire d' Etudes Spatiales et d'Instrumentation en Astrophysique (LESIA). The 
reduction includes the flagging of data obtained during passes of the satellite through the 
South Atlantic Anomaly (SAA), the calculation of a photometric mask, the detection of 
cosmics, contaminants, outliers due to the satellite jitter and hot pixels and filtering of the 
orbital signal. A detailed description will be given in Auvergne et al. (2010, in preparation). 

In preparation for the NGC 2264 observations, a prioritized target catalog consisting 
only of the stars brighter than the Exo-CCD magnitude limit was generated. From the 20 
possible windows on the Exo- field CCD1, 14 were used to observe some of the previously 
selected objects including the two known PMS pulsators V 588 Mon and V 589 Mon. Figure 
Q] shows 5-day subsets of the respective light curves (panels c and f). Note that the different 
filter bandwidths of the CoRoT (from 370 to 1000 nm) compared to the MOST (from 350 
to 750 nm) data are caused by the different passbands used in the two satellites. 



2.3. Frequency Analysis 



Frequency analyses were conducted using Period04 (ILenz fc Breger 1120051 ) which com- 
bines Fourier and least-squares algorithms. Frequencies are prewhitened subsequently and 
are considered to be significant if their amplitudes exceed 4 time s the local noise level in the 



amplitude spectrum (IBreger et al. Ill993l ; iKuschnig et al. 1119971 ). 



For the MOST data sets, the formally significant frequencies were checked against the 
instrumental frequencies related to the orbit of the satellite, its harmonics and 1 d -1 sidelobes 
within the frequency resolution (computed according to iKallinger et al. I l2008bl ). In the 
CoRoT data the influence of the satellite's orbital and related frequencies is negligible. 



3. Pulsation and granulation of V 588 Mon and V 589 Mon 
3.1. Modeling the stellar background in the CoRoT data 



Assuming white, i.e. frequency independent, background noise, the 23.4 days long 
CoRoT time series for V 588 Mon and V 589 Mon showed 106 and 197 formally significant 
frequencies, respectively, after the initial frequency analysis with the methods described 
above. Although the number of frequencies are lower than the many hundreds of f requencies 



reported for classical Scuti stars observed by CoRoT (e.g., iPoretti et al. 1 12009k observed 



>1000 frequencies in HD 50844), they are still too many to be accounted for as simply 
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individually excited low /-valued p-modes. Asteroseismic models predict some tens of p- 
modes with degrees I < 4, which are potentially excited in the frequency range of 5 Scuti- 
type oscillations. A substantially larger number of modes would have to include modes of 
degree I > 4 and/or rotationally split modes and/or the possibility that the detected modes 
are part of the intrinsic background noise due to granulation. 

If all /-valued p-modes are driven to similar amplitudes then the geometric cancelation 
effect reduces the observable amplitude by more tha n an order of magnitude for / > 3 p-mode s 



see, for example, the spatial response fun ctions in IChristensen-Dalsgaard fc Gough Ill982l ). 



But iDaszynska-Daszkiewicz et al. I (120061 ) note, though, that some stars exhibit more than 
two orders of magnitude variations in pulsation amplitudes. In addition, when observing in 
luminosity the geometric cancelation effect is offset by an I 2 factor. Regardless, the variation 
in amplitude of the modes in the short run data (23 d) from CoRoT and in the MOST data 
is a little more than one order of magnitude. We are not yet reaching the lower amplitude 
modes, for which these effects b ecome important. Fu rthermore, we note that unlike the 



example of HD 50844 observed by iPoretti et al. I (120091 ) we do have data for V 588 Mon and 



V 589 Mon from different instruments obtained at different epochs. If all the frequencies were 
due to pulsation modes then they (or at least the majority) should appear in all data sets. 
However, apart from the frequencies listed in Tables [2] and [3] we find 5 to 10 frequencies that 
coincide (within the frequency uncertainties) in the CoRoT and MOST data sets for both 
stars. The number depends on the considered frequency range, the chosen significance limits, 
and how the coincidence and frequency error is defined in detail. But even for two totally 
unrelated data sets it would not be surprising to find a few coincident frequencies. In case of 

V 588 Mon, for example, most of the 106 formally significant frequencies in the CoRoT data 
are distributed between about 5 and 15 d _1 in about 700 independent frequency bins. The 
frequency bin-size is on e over the data set length times the square root of the significance 



( iKallinger et al. Il2008bl ). For our analysis we use an average significance of 10. Even if the 
formally significant frequencies in the CoRoT data set were just random numbers one could 
expect to find by chance matches for at least a few of them with the 179 formally significant 
frequencies determined from the MOST data. Hence, at this time, we argue that many of the 
low amplitude frequencies in the CoRoT data are due to granulation and a (rare) coincidence 
with low amplitude MOST frequencies is not indicative for pulsation. 



Kallinger fc Matthews I (120101 ) concluded that many of the low-amplitude frequencies in 
S Scuti stars are consistent with a strongly frequency- dependent intrinsic background signal 
due to granulation or some effect similar to granulation. This is well known for cool stars with 
convective envelopes. The turbulent motions in their outer convective envelopes generate 
quasi-stochastic power, with amplitudes strongly decreasing for shorter time scales. Although 
the signal is quasi-stochastic in the time domain, its Fourier transform is Lorentzian-like, 
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where different physical processes on (or near) the stellar surface produce the same type of 
signal but on different amplitude and time scales. For the Sun and other Sun-like stars, 
the different signal components are usually assigned to stellar activity, activity of the pho- 
tospheric/chromospheric magnetic network, and granulation, with time scales ranging f rom 
months for active regions to minutes for granulation (see e.g., iMichel et al. 1120081 120091 ). 

The presence of a granulation background signal depends of course on the presence of 
a surface convection zone, which is usually attributed to stars cooler than the red border 
of the classical instability strip. On the other hand, it is known that stellar evolutionary 
models for stars in the instability strip predict a thin convective surface layer and they should 
therefore show a similar signature of granulation as cool stars in their power spectra. We 
believe that this was not recognized before in classical 5 Scuti stars, but also in PMS 5 Scuti 
stars because the amplitudes are small and a clear detection requires long and uninterrupted 
observations to achieve a high-frequency instrumental noise well below lOOppm, which have 
become available only recently. Although the MOST observations of V 588 and V 589 Mon 
are of unprecedented length and completeness the granulation signal with a low-frequency 
amplitude of roughly lOOppm (see Figure [2] and [3]) is hidden in the instrumental noise (see 
Table [T]). Only the CoRoT observations are sufficient to reveal the strongly frequency- 
dependent background signal. 



3.1.1. V 588 Mon 



The CoRoT light curve of V 588 Mon consists of 56978 data points sampled with 32 
seconds and has a noise level of 47 ppm in the original amplitude spectrum from to 100 
d _1 (see Tabled]). The first frequency analysis resulted in 106 formally significant frequencies 
with signal-to- noise v alues larger than 4. Am ong those are also the previously published 8 



pulsation frequencies (IGuenther et al. 1120091 ) 



A visual inspection of the power spectrum of V 588 Mon reveals that the average power 
per frequency bin (e.g. 10 /iHz) is roughly constant at low frequencies and drops by more than 
one magnitude beyond about 500 /iHz, which is a clear indication for a frequency-dependent 
noise. To model the background signal, the pulsation power has to be taken into account as 
well. This is relatively straight forward for cool stars, where the envelope of the pulsation 
power excess is well approximated by a Gaussian. This is more complicated for stars in the 
instability strip where the excess envelope can be quite different from a Gaussia n (or any 
other simple function). We follow the approach of iKallinger fc Matthews I (120101 ) and first 
prewhiten the 10 frequencies with the most significant peaks (i.e., the highest amplitudes) 
and assume the residual pulsation power excess to be at least Gaussian like. Although NGC 
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2264 was observed continuously for the 23.4 days, several data points were removed during 
the reduction process mostly due to high energy particle hits during the satellites passes 
through the South Atlantic Anomaly (SAA) leading to regular gaps in the time series. The 
resulting alias peaks would significantly distort a fit to the power spectrum. We therefore 
filled the gaps - but only for the time series used to investigate the background characteristics 
which requires uninterrupted data - with linearly interpolated values to receive the needed 
clean window function. The power spectrum of the gap-filled data is shown in Figure [2j The 
residual power spectrum is modeled with a superposition of white noise, two Lorentzian- 
like functions, and a Gaussian. The resulting global fit (solid line in Figure [2]) reproduces 
the overall shape of the residual power spectrum and demonstrates the strong frequency 
dependence of the background signal. The residual power spectrum is then corrected for 
the background signal by using the global fit without the pulsation component. Finally, we 
use the height-to-background ratio (HBR; middle panel of Figure [2]) to rate the significance 
of the individual frequencies received from the initial frequency analysis and consider only 
frequencies that exceed an HBR value of 9 to originate from pulsation. For V 588 Mon four 
frequencies match this criterion (see Tableland Figure[2]). Together with the 10 prewhitened 
frequencies of highest amplitudes, a total of 14 pulsation frequencies remain after granulation 
modeling. There are a number of peaks in the middle panel of Figure [2] between about 70 
and 200 /xHz that have a higher HBR value than the peaks in the surrounding frequency 
ranges but which do not exceed our formal limit of HBR > 9. This indicates that there 
are additional pulsation frequencies but our simplistic model of the pulsation power excess 
underestimates their significance. 



3.1.2. V 589 Mon 

The CoRoT light curve of V 589 Mon has 57092 data points also obtained with an 
exposure time of 32 seconds and the noise level in the amplitude spectrum of the original 
data set is 61 ppm from to 100 d _1 (Table [T]). After the initial frequency analysis, 197 
frequencies were identified to be for mally significant. Simil ar as for V 588 Mon, all previously 



published 12 pulsation frequencies ( jGuenther et al. 1120091 ) were also detected in the CoRoT 
data. 

The same method as described for V 588 Mon was used for V 589 Mon. In this case the 
15 frequencies with the highest amplitudes were prewhitened before the gaps in the residuals 
were filled by linear interpolation (Table |3]). The top panel of Figure |3] shows the power 
spectrum of the gap filled data (grey), the global model fit (solid line), the model without the 
pulsational component (dotted line), and the two Lorentzian-like functions (dashed lines). 22 
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frequencies show an HBR value larger than 9 in power (middle panel in Figure [3]). Together 
with the 15 frequencies of highest amplitudes that were prewhitened, a total of 37 frequencies 
are attributed to be caused by pulsation (see Figure Eland Table [3]). 



3.2. Comparison of the MOST 2004 and 2006 data to the CoRoT data 

The MOST 2004 light curves of V 588 Mon and V 589 Mon have 59344 and 60386 data 
points, respectively. The noise levels in the original amplitude spectra computed from to 
100 d -1 (i.e., to 1157 //Hz) are 125 and 133 ppm (see Tabled]). In 2006 V 588 Mon and 
V 589 Mon were observed only less than half the time as for 2004, hence the light curves 
consist only of 23137 and 23113 data points, respectively, and the resulting noise levels from 
to 100 d _1 in the original amplitude spectra are 405 and 429 ppm (Table [1]). 

For each star a separate frequency analysis was conducted and the frequencies common 
to both data sets were then used for the further analysis and for comparison to the CoRoT 
data. The modeling of the background noise was not possible due to the lower quality of the 
MOST data compared to the CoRoT observations. 

For V 588 Mon 20 significant frequencies are common to the MOST data sets of 2004 
and 2006. A cross identification with the 14 CoRoT frequencies that remain after granulation 
modeling showed that 12 of them are also present in the MOST data (Table [2]). A further 
cross-check of the formally significant CoRoT frequencies before the background noise model- 
ing to the 20 frequencies appearing in the MOST 2004 and 2006 data resulted in 7 additional 
common frequencies (F15 to F21 in Table [2]) • Although these peaks have HBR values lower 
than 9, hence would be attributed to be caused by background noise, they have to be con- 
sidered as pulsational. Such stable frequencies that are present in three independent data 
sets obtained in three different years are unlikely to be caused by granulation. Additionally 
the shape of the pulsational power excess i s very likely more complicate d than a Gaussian 



that is used in the power spectrum model (IKallinger &: Matthews Il2010l ) and in the future 



a more realistic approach will be needed. Therefore, for V 588 Mon 19 pulsation frequencies 
are in common between the MOST 2004, MOST 2006 and CoRoT data sets (Table [2]) and 
are marked as black lines in the bottom panel of Figure HJ A comparison between the 106 
formally significant frequencies and the identified 21 pulsation frequencies for V 588 Mon 
is given in the bottom panel of Figure [2] and illustrates the influence of granulation to the 
frequency spectrum. 



For V 589 Mon 19 significant frequencies are common to the MOST data sets of both 
years. 17 of those were also identified in the CoRoT data after granulation modeling and 
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are marked as black lines in the bottom panel of Figure [5j The two additional frequencies 
appearing in the MOST data were identified likely to be aliases or caused by instrumental 
effects. The comparison between the formally significant 197 frequencies and the identified 
37 pulsational peaks is shown in the bottom panel of Figure [31 

Using both MOST data sets and the COROT data, in total 21 pulsation frequencies are 
identified for V 588 Mon (black lines in the top panel of Figure \4ty and 37 for V 589 Mon 
(black lines in the top panel of Figure [3J] where 2 and 20 of them, respectively, are unique 
discoveries obtained with CoRoT. All MOST frequencies have a counterpart in the CoRoT 
observations. 



4. Asteroseismic analysis of V 588 Mon and V 589 Mon 



For any seismic modeling of a star to be successful it is critical that the identified 
frequencies are, in fact, p-modes (or g-modes) intrinsic to the star. Errant or spurious 
frequencies are difficult to isolate by modeling alone. Because viable models are based on 
a large parameter space many possible solutions for various combinations and selections of 
the observed frequencies can be produced. In specific terms, when, as is often the case, it 
is impossible to find a model fit to all the observed frequencies, there usually exist multiple 
possibilities to fit some subset of the observed frequencies and no decisive way to prefer one 
solution over the others. 

CoRoT observations of both V 588 Mon and V 5 89 Mon have yielded a large number of 
frequencies, many of which were also seen by MOST ( iGuenther et al. II2009I ). The two inde- 
pendent observations go a long way toward confirming the intrinsic, i.e., non- instrumental, 
nature of the frequencies. Furthermore, since the MOST and CoRoT observations are sep- 
arated in time by several years, they imply the modes are stable over a period of at least 4 
years. 

Both stars show significant v ■ sin i velocities, high enough, to produce resolvable split- 
tings in the frequency domain. Unfortunately, the mode identification is complicated by 
the fact that the predicted splitting frequencies are comparable to the characteristic (i.e., 
asymptotic large) spacing that separates p-mode frequencies, thereby, making it difficult 
to distinguish the split sidelobes from the parent frequencies. Furthermore, if the stars 
are rapi dly rotating then th e splittings are not linearly spaced about the reference fre 



quencv (|Espinosa et al. 112004 ). Recent studies ( IReese et al. Il2009bt iDeupree fc Beslin 112010 



Ouazzani 1 120081 ) show that non-perturbative multi-dimensional approaches are necessary to 
predict the frequency spacings. Although the direct modeling of these frequencies is compli- 
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cated, iReese et al. I (J2009aJ) show how an empirical formula can be used to fit the frequencies 
assuming the interior rotation profile is simple enough. Needless to say, as long as mode iden- 
tification depends on the interior rotation curve, there will always be some ambiguity in the 
actual identifications. 

Finally, we note that the observed frequencies are low, corresponding to model modes 
near the fundamental frequency of the star. Although the frequencies of these modes are 
more sensitive to the interior structure of the star than higher frequency modes, hence, 
better suited to constraining the evolutionary state of the model, they are also subject to a 
mode bumping-like effect. As a consequence, the regular spacing between adjacent p-modes 
is perturbed, hindering mode identification. 

In the next section we provide a general overview of the modeling methodology and 
theoretical considerations used to interpret the oscillation spectra. In sections 4.2. and 4.3. 
we discuss each star separately beginning with V 589 Mon. 



4.1. Method 



4-1.1. Stellar mode and model calculation 



All of our modeling and seismic analyses are based on a dense set of models that cover 
the HR-diagram. The PMS models start above the birthline on the Hayashi track and extend 
to the ZAMS (zero-age main sequence). Model masses range from 0.81 M to 4.99 M & in 
increments of 0.01 M Q . The models themselve s were constructed using the Yale Stellar 
Evolution Code (YREC: iDemarque et al. 1120081 ) incorporating the latest opacities, nuclear 
physics, and equation of state. Each model is resolved into ~2000 shells with two-thirds 
of the shells covering the outer envelope and Eddington gray atmosphere. We estimate the 
model frequency uncertainty to be 0.1%. The grid density is such that ~2000 models lie 
within 2 -sigma of each star's po sition in the HR-diagram, of the ~400 000 models computed 



in total flGuenther et al. 1120091 ) 



The adiabatic p-mode frequencies f or I = 0, 1, 2, a nd 3 of the models were computed 
using Guenther's stellar pulsation code (IGuenther 112004 ). We have assumed that geometric 
cancelation will occur for higher order p-modes making them more difficult to see above the 
background noise level (see discussion in section 3.1). 
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4-1.2. Mode bumping in PMS stars 

When the oscillation frequencies of our PMS models are plotted in an echelle diagram 
(i.e., frequency modulo the large spacing plotted opposite the frequency of the modes) they 
reveal the vertical alignment characteristic of common /-valued p-modes. But maybe surpris- 
ing at first glance, the lowest radial order, n, modes do not stay aligned with the modes at 
higher frequencies but zigzag. This behavior is common to evolved post-main sequence stars 
and is attributed to the increased density gradients in the region surrounding the isothermal 
helium core near the base of the hydrogen burning shell. T he effect is called mode bumpin g 



and is well documented in texts on stellar pulsation (e.g., ICox I Il980l ; lUnno et al. Ill989l ). 
PMS stars do not show the same large gradients in density because nuclear burning has not 
yet begun. 

Regardless, for some PMS stars there exists a slight peak in the Brunt- Vaisala frequency 
in the deep interior (see Figure for V 589 Mon). This bump is enough to perturb the p- 
mode eigenfunction in the interior and affect its frequency. Note, only nonradial p-modes 
show mode bumping since they need to couple with g-modes (for which I = modes are 
undefined). The perturbation in frequency can be seen in a plot of p-mode frequency versus 
(evolutionary) time as shown in Figure [7] for I = 1 p-modes for models along a 2.65M track 
(corresponding to V 589 Mon, see section 4.2.). The large spacing varies as the star evolves, 
changing inversely with the radius of the star. The n = large spacing deviates slightly 
from the large spacing of the other modes at various times during the evolution of the star. 

Unlike the mode bumping that is seen in post-main sequence stars, the PMS bumping 
does not occur abruptly in time: it is only slowly varying as the model evolves toward the 
ZAMS. Consequently, it is not necessary to increase the resolution of the grid to follow 
this behavior as is the case of post-main sequence bumping. Although it complicates the 
identification of the modes at low frequencies, it also provides an additional feature of the 
oscillation spectrum that can be used to confirm the interior structure and evolutionary state 
of the PMS model. 



4-1.3. Searching for the best fitting model spectra 

To find the be st-fit model oscillation spe ctrum to the observed spectrum we searched 
our grid of models (IGuenther fc Brown 1 120041 ) looking for local minima in \ 2 defined by: 



X 



1 N 

-Y 



j 0" obs,i T 0~ mod,i 
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where v bs,i is the observed frequency for the i th mode, v mo dA is the corresponding model 
frequency, cr & S) j is the observational uncertainty for the i th mode, and iV is the total number 
of modes that match the observed frequencies. We estimate the model uncertainty a mo ^i by 
fitting models to the solar oscillation spectrum ( IGuenther 112004] ) . Here we set it to 0.1% the 
frequency of the mode. 



4-1. 4- Rotation 

As we will describe in sections 4.2. and 4.3., we were unsuccessful in finding p-mode 
frequency fits to all the observed modes. Owing to the not insignificant rotation rates of the 
two stars, we considered the possibility that the oscillation spectrum of each star contains 
rotationally split modes. We included only the lowest order approximations in these attempts 
since full nonlinear rotational splitting computations are too laborious to be applied to a 
grid searching methodology. 

As a first attempt we considered a range of fixed width splittings. These are not model 
constructed splittings that depend on the eigenfunction of the mode and the interior structure 
but are simply a constant frequency added to and subtracted from the model frequency in 
order to simulate, approximately, the expected set of frequency splittings. In other words, 
when searching our grid of models we looked for frequency matches between model and 
observations at v ± mAv, where v is a model computed p-mode, m is the order of the 
splitting, and Av is the splitting frequency itself. 

We al s o com puted rotational splittings for a selected set of models, using the formulation 



of iGough I (1198 ll ). for a solid body rotation curve and for a rotation curve that rises near the 
surface. Here we wanted to see if we could perceive any difference between the two different 
rotation curves, e.g., one curve produced slightly lower y 2 spectrum fits than the other. Note 
that Gough's formula is applicable to slow rotation rates only. 



4.2. V 589 Mon 



4-2.1. Observational Constraints 



V 589 Mon is located just below the birthline in the HR-diagram along our 2.65 M & 
PMS evolutionary track as shown in Figure [H Note that its location in the HR-diagram and 
its association with the cluster NGC 2264 precludes it from being a post-main sequence star. 
The birthline (jStahler I Il983l ; iPalla fc Stahler I Il999l ) corresponding to a mass in-fall rate of 
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~ 10~ 5 M Q /year is shown. Following iGuenther et al. I ( 120091 ) we take the luminosity of V 
589 Mon to be log L/L = 1.58 ± 0.1, the effective te mperature to be 6800 ± 350K, and 
we note that the observed v • sini is 60 ± 10 km s -1 (see 
2, values adopted from 2MASS photometry). 



Kallinger et al. Il2008al their Table 



4.2.2. Mode Analysis 
Figure |9] shows the model computed characteristic spacing (approxi mately equal to the 



average large spacing in the asymptotic limit of large radial order n, iTassoul I Il980l ) as 
function of position in the HR-diagram in the vicinity of V 589 Mon (and V 588 Mon). The 
average small spacing, over the observed frequency range, computed from our model grid is 
1.7 ± 0.3^Hz. 

We plot the observed frequencies in an echelle diagram to reveal the vertical alignment 
characteristic of common /-valued p-modes. Models show that the regular spacing between 
p-modes extends all the way to n — 1 for PMS stars. Based on the model results shown in 
Figure [9] we see that the large spacing for V 589 Mon should lie between 20 /iHz and 26 //Hz. 
The C0R0T frequencies have amplitudes down to 0.27 mmag where the MOST frequencies 
have amplitudes at or above 0.46 mmag. Note that this set of MOST frequenci e s inclu de low 



amplitude frequencies that were excluded in the analysis of IGuenther et al. I (120091 ) where 
only the 8 frequencies common to MOST and ground-based observations were presented and 
analyzed. We proceed with our model analysis using the C0R0T frequencies. 

In an echelle diagram, there is some suggestion for vertically aligned sequences of modes 
in the observations, but the effect is confused by the presence of many additional apparently 
randomly scattered frequencies (see Figure [TO]) . What strikes us most are the many more 
frequencies observed by C0R0T between 50 ^Hz and 150 /zHz than can be accounted for 
from just I — 0, 1, 2, and 3 p-modes. As noted in section 4.1 we believe geometric factors 
and the limited S/N of our data rule out the possibility of identifying higher /-values in our 
C0R0T (and MOST) data. 

We searched through all ~400,000 PMS models in our grid comparing their oscillation 
spectra to the observed spectrum and identified those models that had x 2 < 1- We could 
find no model whose oscillation spectrum fit (with \ 2 < 1) more than 9 of the 37 frequencies 
at a time with / = 0, 1, 2, and 3 modes and which falls within the uncertainty box of V 589 
Mon's position in the HR-diagram. The observed frequencies cannot all be accounted for by 
I = 0, 1, 2, and 3 p-modes alone. 
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4-2.3. Rotation 

To account for the extra modes present in the observed frequency range, we consider 
the possibility that rotationally split frequencies are present. An autocorrelation plot of the 
observed spectrum (Figure [TT1) shows a strong broad peak at 4.1 /iHz (fl), with lesser peaks 
at: 7.9 /iHz (~2fl), 20.2 /iHz (f2 ~ large spacing), 15.8 /iHz (f2-fl), and 24.5 /iHz (f2+fl). 
We speculate that the 4.1 /iHz peak is due to m — ±1 rotational splittings and that the 
7.9 /iHz peak is due to the m = ±2 rotational splittings. Recall that the small spacing, as 
determined from our models is ~1.7 /iHz, hence, cannot explain the 4.1 /iHz peak, unless, 
of course, the deep interior structure of our models is grossly wrong. That the small spacing 
peak is not visible also supports our belief that high /-valued modes are not present in this 
data set. The equatorial velocities, v rot , corresponding to a 4.1 /iHz splitting produced by 
solid body rotation for models near V 589 Mon's location in the HR-diagram are shown 
in Figure [12j The rotation rate implied by the splitting frequency is consistent with the 
observed v ■ sinz of 60 ± 10 km s~ x . 

First to show how rotationally split modes fill up the echelle diagram, we computed 
rotational splittings, assuming solid body rotation (corresponding to ~ 4/iHz splitting), for 
the model loca ted closest to V 589 Mon's HR-diagram position, using the formulation from 



Gough I (jl98ll ). In Figure [TD] we show the CoRoT observations, (a) - (d), along with the 
unsplit model modes, (b) - (d), the m — ±1 splittings, (c) and (d), and the m = ±1, ±2 
splittings. The number density and distribution of model modes are a best match when the 
rotational splittings are included. But, clearly, there is no direct correspondence between 
the model frequencies and the observations. 

We repeated our \ 2 search but used oscillation spectra that include a fixed width rota- 
tional splitting, which we varied from 2.5 /^Hz to 7.5 /iHz in steps of 0.1 /iHz. Restricting 
our search to m = 0, ±1 splittings, we could find, in the best case, matches to 22 of the 37 
frequencies at one time and lie within the uncertainty box of V 589 Mon's position in the 
HR-diagram. The splitting frequency of the best case match is 4.5 /iHz. For m = ±1, ±2 
splittings, we could find model fits to 24 of the frequencies again with a splitting frequency 
equal to 4.5 /iHz. If we assume that we are seeing rotational split modes, the small difference 
between fitting just m = ±1 and m = ±1, ±2 modes suggests that, as expected, m = ±2 
(and m = ±3) rotationally split frequencies are not as prominent in the observations as the 
m = ±1 modes. 

In a crude attempt to show that the model fits do indeed prefer a splitting near 4.5 
/iHz we computed the number of models whose oscillation spectra match exactly 17, 18, and 
19 modes, of the 37 observed, and that also lie inside V 589 Mon's uncertainty box in the 
HR-diagram. When plotted as a function of the splitting frequency, Figure [£31 we see that 
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the number of models that fit the observations peaks for rotational splittings between 4 and 
5 /iHz. The frequency dependence of the number of model fits matches the autocorrelation 
function. There is also a peak at 6.5 /iHz. This peak is likely the combined effect of the I 
= small spacing and the rotational splitting. Although we cannot extract any statistical 
quantities from this plot since the Model Count scale is arbitrary and depends on the density 
of models in our grid, it does show that models that have rotational splittings outside the 
range 4 /iHz to 5 /iHz do not fit the data as well as those within this range of splittings con- 
sistent with our interpretation that the observed spectrum contains sidelobes of rotationally 
split modes. 

That we are not able to fit all of the observed modes with our models may, in part, 
be due to the fact that we assumed the splitting frequency is constant, when, in fact, we 
know that the splitting does depend on model structure, the interior rotation curve, and 
the eigenfunction of the mode itself. To get a sense of how important these effects are we 
computed rotational splittings for a solid body rotation curve and for a rotation curve in 
which the angular velocity is constant throughout the interior then linearly increases by a 
factor a l most a factor of two over the outer 2% radius of the star, using the formula by 



Gough I (119811 ). which is applicable to slow rotation rates. 



In Figure [TH we plot the splitting frequency normalized to the surface rotation rate 
frequency as a function of frequency for I — 1, 2, and 3 p- modes. The rotational splitting 
is constant except at the lowest frequencies where from mode to mode it varies by about 
10%. Coincidentally, this is the frequency range over which the oscillation frequencies for V 
589 Mon are observed. At low frequencies the variation for a 4 /iHz splitting is ~0.4 /iHz. 
Since we have taken our modeling uncertainty to be 0.1%, which corresponds to frequency 
uncertainty of 0.1 /iHz for a 100 /iHz mode, our model fits should be sensitive to the variation 
in splitting with frequency. 

To test if we could see any difference in our \ 2 searches, we computed a new model 
pulsation grid that includes the frequency-dependent splittings for each mode assuming a 
solid body internal rotation curve. For a solid body rotation rate corresponding to a 4.1 /iHz 
rotation rate at the equator of the star we did not find any significantly better or poorer 
model fits than our constant frequency splittings test at 4.1 /iHz. Due to the more laborious 
nature of these more refined model fits and their dependence on the internal rotation curve 
itself, we have not yet pursued more thorough searches over a wide range of rotation rates 
and rotation curves. We believe a different strategy is needed. 
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4.3. V 588 Mon 



4-3.1. Observational Constraints 



V 588 Mon is located just below the birthline in t he HR-diagram a l ong a 2.80 M PMS 
evolutionary track, as shown in Figure EJ Following iGuenther et al. I ( 120091 ) we take the 
luminosity of V 588 Mon to be log L/L = 1.73 ± 0.1, the effective temperature to be T e ff 



7450 ± 350 K, and the observed v • sini is 130 ± 20 km s^ 1 (see iKallinger et al. 112008a . 
their Table 2, values adopted from 2MASS photometry). Except for the higher v ■ sinz the 
star appears to be similar to V 589 Mon. 



4.3.2. Mode Analysis 

Figure [9] shows the expected large spacing as obtained from our model grid in the vicinity 
of V 588 Mon's location in the HR-diagram. Accordingly, we use a folding frequency of 22 
/iHz in our echelle diagrams for V 588 Mon. 

We identify 14 significant frequencies in the C0R0T observations above an amplitude 
of 0.5 mmag. MOST detects 24 significant frequencies above this amplitude. As with the 
M OST data for V 589 Mo n, we have included all significant modes not just those selected 



bv IGuenther et al. I ((2009J) found in both ground and MOST data. The C0R0T and MOST 
sets of frequencies are shown in the echelle diagram of Figure [15j Our model analysis is 
restricted to the C0R0T frequencies. We also plot the I — 0, 1, 2, and 3 p-modes for the 
model in our grid that lies closest to V 588 Mon's location in the HR-diagram. 

There is some suggestion of vertical alignment in the observed frequencies as seen in 
the echelle diagram. But unlike V 589 Mon the distribution and density of modes does not 
point to an obvious over-abundance of frequencies. 

As with V 589 Mon we first sought model fits with I = 0, 1, 2, and 3 p- modes to the 
observed frequencies (excluding rotational splittings). We were only able to find models lying 
with the HR-diagram uncertainty box for V 588 Mon that fit more than 4 (of the 14 C0R0T 
modes) with \ 2 < 1- Indeed, for fits to 5 of the frequencies, no models within 2 sigma of the 
star's HR-diagram location were found. Again we are lead to believe the extra modes are 
rotational splittings. 
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4-3.3. Rotation 

Figure [16] shows the autocorrelation function for V 588 Mon's oscillation spectrum. The 
observed v -sini sets a lower limit on the possible rotational splittings of ~6.7 /aHz (assuming 
that V 588 Mon's location in the HR-diagram is correct) and the breakup velocity (where the 
centripetal acceleration at the equator equals the gravitational acceleration at the surface) 
sets an upper limit on the splitting of ~18 yuHz (again assuming the model parameters based 
on V 588 Mon's location in the HR-diagram are correct). The autocorrelation plot does not 
reveal any obvious splitting peak, although the large spacing frequency does appear to be 
present at ~23 /iHz. The possible splitting peaks are located at 8.3 /iHz, 11.9 /iHz, and 15.1 
/iHz. But clearly none stand out. 

We carried out \ 2 searches with fixed rotational splittings ranging from 5 /iHz to 18 
/iHz but unlike V 589 Mon did not find any preferred frequency splitting or range of split- 
tings. We believe that the rotation rate is, in fact, high enough for nonlinear effects to be 
important. At high rotation rates the centrifuga l and Coriolis force term s become impor- 



tant and introduce non-axisymmetric splittings. [Deupree fc Beslin I (120101 ) have computed 



the oscillation spectra of rapidly rotating, m assive (1OM ) ZAMS stars using a full 2.5D 



evolution code (see also lEspinosa et al. 112004 who first show the complications of rotational 



splittings for rapid rotation). Deupree and Beslin's results are enlightening as they show 
that the autocorrelation function (or, equivalently, as used in their case, a Fourier transform 
of the oscillation spectrum) fails to show rotational splitting frequencies (see their Figure 6) 
already for stars with rotation rates at 10% breakup velocity. The rotational velocity of V 
588 Mon as inferred from its v-sini is minimally at 30% breakup velocity. If we assume their 
results apply to our 2.8M Q PMS star then we should expect to see an autocorre lation plot 



without any predominant splitting peaks (as we do). iDeupree fc Beslin I (120101 ) also note 
that once the rotation rate is above 20% breakup velocity the splittings from different modes 
begin to overlap each other further complicating mode identification. And so, our failure to 
find a model oscillation spectrum, even when including fixed value rotational splittings, may 
simply be a consequence of the far more complex nature of the splittings than our models 
are capable of addressing. 



Rather than attempting to fit the frequencies with the asymptotic-like formula of lReese et al. 



(j2009al ). which they obtained from non-perturbative modeling of frequency splittings in 
rapidly rotating modes for simple rotation curves we are pursuing a different approach that 
does not depend on any prior model or rotation curve assumptions. We will report on our 
results in a future publication. 
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5. Conclusions 

The two pulsating PMS stars V 588 Mon and V 589 Mon were observed by MOST in 
2004 and 2006 for in total ~70 days and in 2008 by the CoRoT satellite for 23.4 days during 
the Short Run SRaOl. Detailed frequency analyses were conducted for all available data sets 
and compared to each other. 

Our analysis illustrates that the frequency dependent intrinsic background signal, i.e. 
granulation, can explain a large number of significant peaks detected in PMS stars. Granu- 
lation modeling was conducted using a first order model where two Lorentzian-like functions, 
white noise and a Gaussian are combined. The resulting number of frequencies is more con- 
sistent with the expected number of low-degree p-modes observed in integrated light. But 
it was also shown that the shape of the pulsational power excess is very likely to be more 
complicated than a Gaussian and that a more sophisticated approach would be needed in the 
future. This effect is illustrated by the fact that for V 588 Mon 7 frequencies at low ampli- 
tudes would have been suppressed by the background noise model, but appear significantly 
in all three satellite data sets, i.e. in data from MOST 2004, MOST 2006 and CoRoT. If 
frequencies are stable over a period of 4 years, they are unlikely to be caused by granulation. 

After a comparison of the independent analyses of the MOST 2004 and MOST 2006 data 
sets to the CoRoT data, for V 588 Mon and V 589 Mon 21 and 37 frequencies, respectively, 
can be attributed to p ulsation, among those the 8 and 12 previously published frequencies 



flGuenther et al. 1120091 ) 



Even after granulation filtering, both V 588 Mon and V 589 Mon have more frequen- 
cies observed than can be accounted for by I — 1, 2, and 3 p-modes. We investigated the 
possibility that the extra frequencies are rotational split modes. For V 589 Mon the auto- 
correlation of the observed frequency spectrum and our model searches support the notion 
that rotational split frequencies, with splittings between 4 /zHz and 5 /zHz, are present. We 
are able to match more than half of the observed frequencies with model fits. For V 588 
Mon the autocorrelation plot and our attempts to fit model spectra that included splittings 
were unsuccessful. We believe this is because V 588 Mon's rotation rate is high enough 
as implied by i t s v ■ smi) that its splittings are in the nonlin e ar regime as describ ed by 



Deupree fc Beslin I fcOlOf ) 



Espinosa et al. I (|2004[). a nd further studied by lOuazzani I (120081 ). iReese et al. I ( I2009bl ). and 



How can we obtain unambiguous mode identifications for V 588 Mon and V 589 Mon 
when we do not know the interior rotation curve? Is inversion possible, i.e., determining the 
rotation curve from the splittings, simultaneously with identifying the modes? Unlike the 
slowly rotating Sun, where the splittings are in the linear regime, V 588 Mon and V 589 Mon 
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are rotating rapidly enough that non-perturbative models are probably necessary, especially 
for V 588 Mon. Rather than to continue with more sophisticated models we are currently 
pursuing a different approach in which we apply a Bayesean approach that incorporates 
model independent prior knowledge about the split and unsplit modes to estimate model fit 
likelihoods. We will report on our efforts in a future paper. 
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Fig. 1.— 5-day subsets of the MOST 2004 (a,d), MOST 2006 (b,e) and CoRoT (c,f) light 
curves for V 588 Mon (left) and V 589 Mon (right) to the same y-axis scale. The different 
amplitudes in the MOST versus the CoRoT light curves is mainly due to the different 
bandpasses used in the two satellites. 
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Table 1. Properties of the V 588 Mon and V 589 Mon data sets obtained from MOST in 

2004 and 2006 and from COROT in 2008: time base (T), frequency resolution (fres), 
number of data points (N), point-to-point scatter (pt2pt), standard deviation (sigma) and 
noise level computed from the original data sets from to 100 d _1 (noise level). 



data set 


star 


T 


fres 


N 


pt2pt 


sigma 


noise level 






[d] 


[d- 1 ] 


# 


[mag] 


[mag] 


[ppm] 


MOST 04 


V 588 Mon 


48.13 


0.02 


59344 


0.008 


0.015 


125 




V 589 Mon 


48.13 


0.02 


60386 


0.009 


0.018 


133 


MOST 06 


V 588 Mon 


22.77 


0.04 


23137 


0.013 


0.022 


405 




V 589 Mon 


22.77 


0.04 


23113 


0.013 


0.024 


429 


CoRoT 


V 588 Mon 


23.41 


0.04 


56978 


0.003 


0.007 


47 




V 589 Mon 


23.41 


0.04 


57092 


0.0009 


0.010 


61 
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Fig. 2. — Top panel: Residual gap-filled power spectrum of V 588 Mon (grey) after prewhiten- 
ing the 10 highest amplitude frequencies. The global fit is shown as solid line, while the 
dotted line illustrates the model without the pulsational components consisting of the two 
Lorentzian-like functions and white noise, and the dashed lines mark the two Lorentzian-like 
functions. Middle panel: Residual power spectra normalised to the background components, 
where the dashed line marks a height-to-background ratio (HBR) of 9 (given in power). Bot- 
tom panel: Comparison between the 106 formally significant frequencies (positive values) to 
the identified 21 pulsation frequencies (negative values) which consist of the 10 prewhitened 
peaks before granulation modeling, 4 peaks exceeding an HBR value of 9 and 7 peaks com- 
mon to all available data sets. 
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Fig. 3. — Top panel: Residual gap-filled power spectrum of V 589 Mon (grey) after prewhiten- 
ing the 15 highest amplitude frequencies. The global fit is shown as solid line, while the 
dotted line illustrates the models without the pulsational component consisting of the two 
Lorentzian-like functions (dashed lines) and white noise, and the dashed lines mark the two 
Lorentzian-like functions. Middle panel: Residual power spectra normalised to the back- 
ground components, where the dashed line marks a height-to-background ratio (HBR) of 
9 (given in power). Bottom panel: Comparison between the 197 formally significant fre- 
quencies (positive values) to the identified 37 pulsation frequencies (negative values) which 
consist of the 15 prewhitened peaks before granulation modeling and the 12 additional peaks 
exceeding an HBR value of 9. 
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Table 2. V 588 Mon: 21 pulsation frequencies (freq) identified from the CoRoT data with 

their respective last digit errors given in parentheses, amplitudes (amp), 
height-to-background values (HBR) where appropriate, flag for frequencies also present in 
the MOST 2004 and 2006 data s ets (MOST), and the c ross-identification to the previously 
published 8 frequencies (Gue09; iGuenther et al. 1120091 ). The first 10 frequencies with the 
respective highest amplitudes were prewhitened before granulation modelling. 



No. 


freq 


freq 


amp 


HBR 


MOST 


Gue09 


# 


[d- 1 ] 


[/xHz] 


[mmag] 


power 


y / n 


# 


Fl 


5.1394(7) 


59.484(8) 


5.804 




y 


Fl 


F2 


9.1579(9) 


105.99(1) 


3.853 




y 


- 


F3 


5.2628(8) 


60.91(1) 


3.653 




y 


F3 


F4 


7.029(1) 


81.35(1) 


2.315 




y 


F4 


F5 


10.935(1) 


126.56(2) 


1.779 




y 


F2 


F6 


5.546(1) 


64.19(2) 


1.673 




y 


F5 


F7 


6.213(1) 


71.91(2) 


1.624 




y 




F8 


8.290(2) 


95.95(2) 


1.446 




y 


F6 


F9 


12.930(2) 


149.65(2) 


1.279 




y 




F10 


7.668(2) 


88.74(2) 


1.152 




y 


F7 


Fll 


8.438(2) 


97.67(2) 


1.046 


17.6 


y 




F12 


11.161(3) 


129.18(4) 


0.670 


12.6 


n 




F13 


8.145(3) 


94.27(4) 


0.610 


11.0 


y 




F14 


10.400(4) 


120.37(4) 


0.571 


10.6 


n 




F15 


9.672(4) 


111.95(5) 


0.464 


7.6 


y 




F16 


10.198(6) 


118.03(7) 


0.344 


6.1 


y 




F17 


10.981(6) 


127.10(7) 


0.317 


6.7 


y 


F8 


F18 


12.128(7) 


140.37(8) 


0.280 


5.1 


y 




F19 


5.115(8) 


59.20(9) 


0.243 


4.3 


y 




F20 


11.661(9) 


134.9(1) 


0.203 


4.1 


y 




F21 


15.36(1) 


177.7(1) 


0.175 


4.4 


y 
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Table 3. V 589 Mon: 37 pulsation frequencies (freq) identified from the CoRoT data with 

their respective last digit errors given in parentheses, amplitudes (amp), 
height-to-background values (HBR) where appropriate, a flag for frequencies also present in 
the MOST 2004 and 2006 data se ts (MOST) and the cr oss-identification to the previously 
published 12 frequencies (Gue09: iGuenther et al. 1120091 ) The first 15 frequencies with the 
respective highest amplitudes were prewhitened before granulation modelling. 



No. 


freq 


freq 


amp 


HBR 


MUb 1 


Gue09 


# 


rj — li 


r, u 1 


[mmag] 


power 


y / n 


# 


Fl 


6.4884(5) 


75.097(6) 


11.615 




y 


Fl 


F2 


6.9862(6) 


80.859(6) 


7.551 




V 

J 


F2 


P3 


7.5426(8) 


87.299(9) 


3.640 




n 




Fl 


8.3005(9) 


96.07(1) 


2.954 




v 

J 


F3 


F5 


5.7643(9) 


66.72(1) 


2.719 




v 

J 


Fl 


F6 


8.698(1) 


100.68(1) 


2.048 




V 

J 


F5 


F7 


10.192(1) 


117.97(1) 


1.773 




V 
J 


F6 


F8 


9.467(1) 


109.57(1) 


1.720 




V 
J 


F7 


F9 


9.314(1) 


107.80(2) 


1.299 




V 
J 


F9 


F10 


6.281(1) 


72.70(1) 


1.287 




V 


F8 


Fll 


7.407(1) 


85.73(2) 


1.147 




11 




F12 


6.699(1) 


77.54(2) 


1.063 




V 

J 


F12 


F13 


8.490(1) 


98.26(2) 


1.036 




11 




F14 


10.809(1) 


125.10(2) 


1.007 




11 


- 


F15 


10.990(1) 


127.20(2) 


0.945 




y 


Fll 


F16 


11.043(1) 


127.81(2) 


0.816 


25.6 


n 




F17 


11.372(2) 


131.62(2) 


0.706 


21.6 


y 




F18 


7.876(2) 


91.16(2) 


0.696 


20.5 


y 




F19 


6.138(2) 


71.04(2) 


0.670 


19.3 


y 




F20 


10.697(2) 


123.81(2) 


0.543 


17.0 


n 




F21 


7.205(2) 


83.40(2) 


0.540 


16.1 


n 




F22 


10.262(2) 


118.78(2) 


0.517 


15.7 


y 


F10 


F23 


8.386(2) 


97.06(2) 


0.501 


15.2 


n 




F24 


11.119(2) 


128.70(2) 


0.468 


13.8 


n 




F25 


8.953(2) 


103.62(2) 


0.456 


13.4 


n 




F26 


7.503(2) 


86.84(2) 


0.455 


13.6 


n 




F27 


12.980(2) 


150.23(2) 


0.442 


15.3 


n 




F28 


9.698(2) 


112.24(2) 


0.437 


13.5 


y 




F29 


9.260(2) 


107.17(2) 


0.431 


12.8 


n 




F30 


10.109(2) 


117.00(2) 


0.408 


12.6 


n 




F31 


6.524(2) 


75.51(2) 


0.405 


11.4 


n 




F32 


7.579(2) 


87.72(3) 


0.333 


10.0 


y 




F33 


12.257(2) 


141.87(3) 


0.314 


10.2 


n 




F34 


7.313(3) 


84.65(3) 


0.313 


9.1 


n 




F35 


10.643(3) 


123.18(3) 


0.290 


10.1 


n 




F36 


6.199(3) 


71.74(3) 


0.270 


11.1 


n 




F37 


13.148(3) 


152.18(3) 


0.267 


9.5 


n 
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Fig. 4. — Comparison of the V 588 Mon amplitude spectra (grey) derived from CoRoT 
(oriented upwards) and MOST 2004 (negative values assigned) data where the respective 21 
and 19 pulsation frequencies are identified in black; the dashed line marks the MOST orbital 
frequency and the dotted lines are the respective 1 d _1 sidelobes. 
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Fig. 5. — Comparison of the V 589 Mon amplitude spectra (grey) derived from CoRoT 
(oriented upwards) and MOST 2004 (negative values assigned) data where the respective 37 
and 17 pulsation frequencies are identified in black; the dashed line marks the MOST orbital 
frequency and the dotted lines are the respective 1 d _1 sidelobes. 
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Fig. 6. — Propagation diagram for the interior model of V 589 Mon. The square of the Brunt 
Vaisala and Lamb frequencies, in units of radians 2 s~ 2 , are plotted against radius fraction. 
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Fig. 7. — The time evolution of the large frequency spacing of the lowest n-valued, / = Op- 
modes for a 2.65M PMS model track (which passes through V 589 Mon's position in the 
HR-diagram). The frequency difference between the n — 1 and n = p- modes is drawn 
with a solid line (indicated in the legend by v\ — uO). The frequency differences between 
the n = 2 and n — 1, and n = 3 and n = 2 are drawn with dashed lines (as labeled in the 
legend). The birthline, the ZAMS, and the point where the track passes through V 589 Mon 
are indicated. 
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Fig. 9. — Contour plot of characteristic spacing (asymptotic large spacing), in units of /iHz, 
computed from PMS models near V 588 Mon's and V 589 Mon's locations in the HR-diagram 
as indicated by the uncertainty bars. 
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Fig. 10. — Echelle diagrams for V 589 Mon with folding frequency 20.3/iHz. Each panel adds 
in sequence, starting with the CoRoT observed frequencies in (a), the model I = 0, 1, 2, and 
3 p- modes in (b), the m = ±1 model rotational split frequencies in (c), and the m = ±2 
model rotational split frequencies in (d). 
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Fig. 11. — Autocorrelation plot of the oscillation frequency spectrum of V 589 Mon observed 
by CoRoT with the frequencies of the highest peaks labeled. 
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Fig. 12. — Equatorial rotational velocities in units of km s 1 are shown in the vicinity of V 
589 Mon's HR-diagram location. The rotational velocities are computed (they depend on 
the model radius) for a rotational splitting (solid body) of 4.1 /iHz. 
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Fig. 13. — The autocorrelation plot for CoRoT observations of V 589 Mon (dashed line) 
compared to the number of models that fit 17, 18, and 19 of the 37 observed frequencies for 
a given autocorrelation frequency (//Hz), and that lie within the uncertainty box of V 589 
Mon's position in the HR-diagram. The lowest curve (solid line joining points) corresponds 
to the number of models fitting 19 frequencies. The next lowest curve corresponds to the 
number of models fitting 18 frequencies and the highest curve to the number of models fitting 
17 frequencies. 
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Fig. 14. — Normalized rotational splittings (Au/u) as a function of frequency in /iHz, v, 
computed from a model of V 589 Mon for solid body rotation (solid line) and for the rotation 
curve shown in the insert (dashed line). 
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Fig. 15. — Echelle diagram showing all the significant frequencies observed by CoRoT and 
MOST of V 588 Mon falling between 50 ^Hz and 175 ^Hz. The folding frequency is 22//Hz. 
The / = 0, 1, 2, and 3 p- modes of the model that falls closest to V 588 Mon's position in 
the HR-diagram are indicated by a solid line connecting large points, a solid line connect- 
ing smaller points, a short dashed line connecting smaller points, and a long dashed line 
connecting smaller points, respectively. 
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Fig. 16. — Autocorrelation plot of the oscillation frequency spectrum of V 588 Mon observed 
by CoRoT with the frequencies of the highest peaks labeled. The lower limit to the splitting 
frequency set by the observed v • sini and the splitting frequency corresponding to the 
equatorial breakup velocity are also indicated. 



